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ABSTRACT 


In an earlier thesis, Ellsworth showed that the radiated power output (quantified by 
the gain in radiation resistance presented to the transducer) and directivity of a compact 
underwater transducer can be significantly improved by the use of a system of resonant 
scatterers. These resonant scatterers were termed "sympathetic resonators." In the 
present work, we verify the previous findings and extend the work to additional cases of 
interest. Specifically, we (1) employ equations requiring fewer assumptions and 
implement them using more accurate numerical techniques, (2) reproduce Ellsworth’s 
calculations and provide graphical results for conical and planar circular configurations 
for the case of six passive acoustic resonators, (3) extend the theory to and provide 
results for two types of linear arrays of passive acoustic resonators, (4) and finally, we 
consider an acoustic Yagi-Uda antenna-like configuration of passive acoustic resonators 
to improve the directivity of an underwater transducer. The results indicate that the 
conical configuration of six resonators is the best choice in terms of the gain in radiated 
power output, directivity along an acoustic axis, and minimum required number of 


resonators. 
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I. INTRODUCTION 


A. BACKGROUND 

It has been shown that a coupled set of resonant scatterers will, at frequencies near 
the resonance of a single scatterer, exhibit a quasiresonant condition, resulting in a 
significantly higher scattering cross section than a single scatterer [Ref. 1], [Ref. 2]. 
This phenomena has been shown by I. Tolstoy to provide significant gain in the scattered 
acoustic pressure to a plane wave when incident upon a system of scatterers as in Figure 


1-1 [Ref. 3]. Tolstoy demonstrated that, particularly for the line array, the amplification 


Plane Wave 





Resonant Scatterers 


Figure 1-1 Plane wave incident upon a system of resonant scatterers. 


of the scattered wave is sensitive to the direction of the incoming plane wave (Ref. 3]. 
In his thesis [Ref. 4], J. Ellsworth considered the effect of the presence of a system of 


resonant acoustic scatterers on a transducer in the coupled acoustic field as in Figure 1-2. 


Ellsworth referred to these resonant scatterers as "Sympathetic resonators." We will 
refer to these scatterers as passive acoustic resonators. Two geometries were previously 
considered. The first, called the planar configuration, placed the transducer in the center 


of a circular array of resonators. The second, called the conical configuration, is the 


Resonant Scatterers 
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Transducer 


Figure 1-2 A transducer in a coupled field of resonators. 


same as the planar configuration except that the transducer is displaced one-quarter 
wavelength out of the plane of the resonators. The results demonstrated a gain in 
radiation resistance presented to the transducer. This means that, for a displacement 
limited transducer, more power can be radiated into the medium. Ellsworth also 
provided beam patterns demonstrating that the transducer will, in the presence of passive 


acoustic resonators, exhibit directionality along an acoustic axis. 


B. OBJECTIVES 
The objectives of this thesis are to (1) refine the theory presented by Ellsworth 
using equations requiring fewer assumptions and implementing those equations with more 


accurate numerical methods, (2) reproduce Ellworth’s work using the more accurate 


techniques and provide three-dimensional plots for the case of six resonators (the six 
resonator results were requested by the Naval Ocean Systems Center), (3) develop three 
dimensional beam patterns for all cases of interest, (4) develop the theory for and provide 
results for two types of line arrays of passive acoustic resonators, and (5) to consider an 
acoustic Yagi-Uda [Ref. 5] antenna-like configuration of scatterers to improve the 
directivity of an underwater transducer. 

Our measures of effectiveness in the above work will be (1) the gain in radiation 
resistance presented to the transducer with the system of passive acoustic radiators 
present, (2) three-dimensional plots of that gain versus the system parameters, (3) the 
directivity index of the transducer in the presence of passive acoustic radiators, and (4) 
the shape of and characteristics of the far-field beam pattern for each configuration 
considered. In support of these objectives, programs were written in Mathematica [Ref. 
6] (a system for doing mathematics oy computer) on Macintosh, IBM, and SUN 
workstation computers. An example program is provided in the appendix. 

The remainder of this thesis is organized as follows. We first develop the theory 
for a transducer in a coupled acoustic field with passive acoustic resonators. We provide 
results for the two circular configurations, planar and conical. We analyze and provide 
results for two types of line arrays of resonators, termed line type I and line type II. We 
consider the acoustic Yagi-Uda array of resonators. Finally, we discuss the results and 


provide conclusions and recommendations. 


Il. THEORY 


A. THE NETWORK MODEL 

We will develop the network model for a coupled network of a transducer and N 
passive acoustic resonators by considering each element as a lumped acoustic system. 
The development is similar to that presented by Ellsworth [Ref. 4]. We will represent 
each lumped acoustic element as an acoustic impedance, choosing pressure and volume 
velocity as mechanical variables. In Section 1 we develop the network representation of 
a single resonator in an acoustic field. In Section 2 we develop the network 
representation of a transducer in an acoustic field. In Section 3 we combine the results 
of Sections 1 and 2, and develop a network model for a transducer coupled to N 
resonators by an acoustic field. We allow for interaciions between each element and 
every other element. This network will be used to analyze systems of passive acoustic 
resonators in the chapters that follow. Finally, in Section 4, we discuss the required 


input parameters for the model developed. 
1. A RESONATOR IN AN ACOUSTIC FIELD 
Figure 2-1 shows a lumped acoustic model of a single resonator in an acoustic 
field. The network equations for this system are: 
p, = -Zf U, (2-1а) 


апа 


Р! x Zn U, рү (2-10) 


where p, = acoustic pressure at the surface of the resonator 


- 


Resonator Acoustic Field 





Figure 2-1 Lumped acoustic model of a resonator in an acoustic field. 


U, = acoustic volume velocity of the resonator 
Z° = acoustic impedance of the resonator 
z% = acoustic radiation impedance presented to the resonator by the 


medium with all other active or passive elements blocked 


рр = the acoustic pressure acting on the surface of the resonator when 


it is blocked (1.e. U, = 0) 
D = diffraction constant [Ref. 7] 


Du = free-field pressure incident on the resonator, i.e. the pressure that 


would exist at the location of the resonator if 1t were not present. 


where 


The free-field pressure can be found from 


рї = Le zu, (2-2) 
m increments over all other active and passive elements present and 
Zi” = the free-field acoustic transfer impedance from element m to the 
resonator 
U, = acoustic volume velocity of element m. 


Bobber [Ref. 7] defines the diffraction constant as the following ratio 


р. = P (2-3) 


The blocked pressure is the average pressure over the surface of the resonator if its 


volume velocity is zero. The surface pressure 1s, in general, the sum of the blockea 


pressure and the radiated pressure, as indicated in Figure 2-1. The radiated pressure is 


given by 


рі“ = iy (2-4) 


and so the total pressure at the surface is 


p, = pi +p“ =Dpf+ UZÎ . (2-5) 


For notational convenience, we will include the diffraction constant as part 


of the acoustic transfer impedance between two elements. Let Z^ be the transfer 


impedance from element m to element n. We can then write 
ees L zuy. (2-6) 


The radiation impedance is a function of the properties of the medium and 
the physical size and shape of the object radiating into the medium. The input acoustic 


impedance is a function of the properties of the resonator. 


2. A TRANSDUCER IN AN ACOUSTIC FIELD 


Figure 2-2 is the lumped acoustic model of a transducer in an acoustic field. 


Transducer Acoustic Field 














Figure 2-2 Lumped acoustic model of a transducer in an acoustic field. 


The network equations for this system are 


Die Тыс де U, (2-7а) 


Po = Zoo Ue Po (2-7b) 


and 


where p. 
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ла 


lo 


V, 


| Z 


EB 


1 


We E (2-7с) 
acoustic pressure at the surface of the transducer 


acoustic volume velocity of the transducer 
open-circuit acoustic impedance of the transducer 


transduction coefficient of the transducer 
acoustic radiation impedance presented to the resonator by the 


medium with all other active or passive elements blocked 


the blocked acoustic pressure at the surface of the transducer 


blocked electrical impedance of transducer 
voltage across the electrical terminals of the transducer 


current through the electrical terminals of the transducer. 


The blocked acoustic pressure is found as before. The electrical quantities 


and physical characteristics of the transducer will not be considered further here but are 


included for completeness. The transducer is taken to be compact and we are interested 


in the radiation impedance presented to it. The volume velocity of the transducer will 


be held constant (assumed to be displacement-limited). For a displacement-limited 


transducer, the power radiated can only be increased by increasing the radiation loading. 


3. A TRANSDUCER IN THE PRESENCE OF PASSIVE ELEMENTS 
Figure 2-3 is the network model of N interacting passive acoustic resonators 


coupled to a single transducer by an acoustic field. This model combines the results of 
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Figure 2-3 Network model of N interacting passive acoustic resonators and a transducer 


the previous two sections and allows for interactions between each element and every 


other element. The network equations are 


IS LA (2-8a) 


and 


Pr 2 n=0,1,2,.....N (2-8b) 


p, = 2 2 NENN (2-8с) 
where N = number of passive acoustic resonators 
p, = acoustic pressure at the surface of element n 
U, = volume velocity of element n 
Z7 = acoustic radiation impedance seen by element n with all other active 
or passive elements blocked 
Z. = acoustic transfer impedance from element m to element n. 


We are interested in the total acoustic radiation impedance seen by the 


transducer in the presence of the passive acoustic resonators. This quantity can be 


defined as follows 


| (2-9) 
9 


Equation 2-8 can be written for the case of N=6 resonators as the following 


matrix equation 


10 


ar ar ar ar ar ar ar 
Lo Zo Lo Zo Zo Zo Los 


Po Uo 

E Z5 Zi Zu Za Za Z5 Z«| |у 

р, Us Ui Ue Ue тк шш 
EH Т Za ZE AU, SE 

Р, ZA. 74877. ZO n Ox 

ar ar ar ar ar ar ағ U 

Ps WC ZU Zum Zu s 

Ps | ar ar ar ar ar ar ar U, 

“о Zei Lo Lo Ze Los Los 
In order to simplify the notation, we will write equation 2-10 as 

Po - Zoo 24] D: (О 1) 


|Р.) Eh kl QU 


where x) is a IxN column matrix 


[x] is a NxN square matrix 
and the superscript "7" indicates the transpose of a matrix. 


In order to find Z7, first we replace {р by |-ZzU,| and normalize by U,, 


total ? 


resulting in the following equation: 


ar ar ar 7 
Ze] | Za (е) С | us 
(0) (zh Ie +z? oil Lom 


where б, = | 
апа US ae ae 


We can represent the matrix of equation 2-12 as the following system of equations: 


2 T NES Z luz (2-13) 
and . 
lo) = {25} + [25 + 722 UY M 


Solving por) in equation 2-14 yields 


ful = -[ze+z2 5) bh (2-15) 


Equation 2-15 gives the relative volume velocities of each resonator with respect to the 


transducers volume velocity. We substitute equation 2-15 into equation 2-13 to obtain 


the following expression for 24”; 


zi, = Zu - bal el bel, (2-16) 


The total radiation impedance presented to the transducer will be used to find the gain 


in radiation resistance with resonators present. 


4. THE INPUT PARAMETERS 
The physical parameters that we must provide in the model developed above 
are (1) the acoustic impedance of each resonator, (2) the transfer impedance of each 
element in the presence of other elements and with the other elements blocked, (3) the 
physical characteristics of the medium (hydrostatic pressure, sound speed, etc.). The 
required equations are developed later in this chapter. The equations for the measures 
of effectiveness chosen to evaluate the resonator systems are developed in the sections 


that follow. 


B. GAIN IN RADIATION RESISTANCE 


The total power radiated into the water by a simple source is [Ref. 8] 


ne — | * Re[f] Re[ugldt (2-17) 
Ша 
where Da force of the fluid on the source 


particle velocity at the surface of the source 


A 
о 
| 


ل“ 
| 


period of harmonic motion. 


In terms of the variables used in this paper, the total power delivered to the acoustic field 


by a transducer in the presence of N passive resonators is given by 


1 аг 
П, = А OREZ (2-18) 


The gain in power delivered to the acoustic field due to the presence of the 


resonators can be defined as [Ref. 4] 


D 


total power delivered with the resonators present (2-19) 


Gain LLLI eee 
total power delivered without the resonators 


For a compact source, the radiated acoustic power in the absence of the resonators can 


be approximated by 


IL = 5 08 25). (2-20) 
Then we can write equation 2-19 as 
Kess 
Саїп = Кешш (2-21) 
Ке[20] 


We will refer to this quantity as the gain in radiation resistance, or simply gain. The 
gain in radiation resistance 1s an important measure of effectiveness because it quantifies 
the additional radiated power due to the presence of passive acoustic radiators. We will 
use this quantity to evaluate the benefit of using passive acoustic radiators in various 


configurations. 


C. BEAM PATTERNS AND DIRECTIVITY 

The directivity of a transducer due to the presence of passive acoustic resonators 
can be determined by analyzing the far-field beam pattern. To find the beam pattern we 
find the radiated pressure due to each element at far-field points in various directions. 


For a compact source, we can calculate the field pressure due to a single element as 


14 


-jkr 
а, е JH, 


б слу к. = (2-22) 
г. өзін, 
where pl" = acoustic pressure at a field point due to element n 
a, = radius of element n 
r = distance to the field point from the acoustic center of element n. 


The total pressure at the field point due to the transducer and system of resonators is then 
N 
pa (r8,6) болад) (2-23) 
n=0 


where р“ (т,@,фу = the pressure at the field point specified by the spherical 


coordinates (r,9,4) due to element n. 
In the plots which will be presented, the above far-field pressure is normalized by 
the pressure at the same point due to the transducer alone. What is plotted, then is the 


following function, termed BP(@,¢): 


im | pe 6.9 | (2-24) 


BP(@,¢) = 
IN МЕС | 


We will plot the beam pattern defined by equation 2-24 in three dimensions in order to 
easily see what is the effect of the resonators on the pressure field. We will also plot the 


beam pattern on a log scale in the form of a polar plot. Polar plots are the traditional 


15 


method of presenting a beam pattern. In this case 0 dB will represent the sound pressure 
level of the transducer with no resonators present. 
A measure of how well the system projects power along an acoustic axis is the 


directivity index. The directivity index is defined as [Ref. 8] 








DI = 10 oe Бе (2-25) 
p; (7) 
where р.) = acoustic pressure due to the directional source at a range r 
along the acoustic axis 
joue cm acoustic pressure, at a range r, due to a spherical source 


radiating the same acoustic power as the directional source. 


In terms of our equations, the directivity index can be show to be 
Di =" 20 Jog (ВИО а 20108 (Gain) Gm 


where 0, and ¢,, indicate the spherical angles along the acoustic axis. 


D. GAS-FILLED BUBBLES AS PASSIVE ACOUSTIC RADIATORS 
Tolstoy and Ellsworth both considered air-filled bubbles as passive acoustic 


resonators in their papers. Here, we will consider air-filled bubbles at a hydrostatic 


pressure Py — 1 atm (just below the surface). We will assume a sound speed in water 


16 


of c = 1500 meters per second and a density in water of p, = 998 kg/m’. The ratio of 
specific heats for airis y = 1.4. 

In addition to the above constants, we require the acoustic impedance of each 
resonator and the transfer impedances between all elements. The required impedances 
for air-filled bubbles are developed below. 


The mechanical stiffness of an air bubble in water is [Ref. 8] 
рКа т (2-27) 


where a is the radius of the bubble. This stiffness can be represented as an acoustic 


impedance as follows 


eo (2-28) 
S ca 
where $ = surface area of the bubble 
k = wavenumber in water. 


We have only included the adiabatic stiffness of the air, since the effects of surface 
tension, non-adiabatic expansion and contractions, and losses due to viscosity can be 
neglected for bubbles larger than 1 mm [Ref. 9] and frequencies below 40 kHz [Ref. 8]. 


The resonant frequency of a bubble in water is [Ref. 8] 


ЖҮР, (2-29) 


GQ = 1 
0 
a Po 





We will prefer to consider dimensionless quantities such as the wavenumber times 


distance. Equation 2-29 can then be expressed as 


(2-30) 





For the values specified above, the resonant value of ka is ka, =0.01377. 


The acoustic self-impedance seen by element n is [Ref. 8] 


poo ka, 
5 y1 + (ka, 


where cot(0) = ka,. This can be easily rewritten as 


p. Өй (2-31) 


zz = Ро Ха, ` (2-32) 
S 1+jka, 





The transfer impedance from element m to element 7 1s 


ze = Ze D "pits, - I) (2-33) 


Kl. 
where e: is the distance between the acoustic centers of elements m and n. The 
diffraction constant for a sphere of radius a, is (Ref. 7] 

jka, 


= ШШ (2-34) 
1 +jka, 


л 


We have taken the transducer radius, dy, to satisfy kag = 0.2. This is a typical 
value for a low frequency transducer. Later we will show that the results are not 
sensitive to this choice over a large range of values. This value of transducer radius was 


used for all calculations unless otherwise noted. 


E. OTHER PASSIVE ACOUSTIC RADIATORS 

The network model developed above is general and can be used for a variety of 
passive acoustic radiators. Any compact resonator that can be characterized by a 
stiffness alone, such as an air bubble in water, can use the impedances developed for the 
bubble with the appropriate choice of ka,. In addition to the air bubble at 1 atm ambient 
pressure, we have chosen to perform calculations for a bubble-like resonator with a value 
of ka, = 0.1 (which corresponds to an air bubble at 50 atm ambient pressure). We can 
expect that a very thin-shelled object made of steel or aluminum, an air balloon, or air 
bladder might be built with these characteristics. We have not attempted to devise or 
construct such a system, but rather present results to determine whether one should be 
built. We will compare the results for air bubbles at 1 atm ambient pressure and 


resonators with ka, — 0.1. 
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III. THE CIRCULAR CONFIGURATIONS 


In this chapter, we present the results for the two circular configurations previously 
analyzed by Ellsworth [Ref. 4]. The tables of results presented below were created using 
more accurate numerical techniques than those used by Ellsworth. Also, the equations 
used require fewer assumptions. Specifically, (1) we use an exact equation for the 
radiation impedance from a sphere (rather than an approximate formula), (2) we take into 
account the diffraction constant relating the blocked pressure to the incident pressure 
wave, (3) we used a gradient search method to find the maximum gain as a function of 
the free parameters rather than reporting the maximum found over the coarse grid of 
parameter values, (4) we calculate the directivity index for each case considered, and (5) 
we plot three-dimensional beam patterns in order to better visualize the effect of passive 
acoustic resonators on the acoustic field. Ellsworth provided plots of gain in radiation 
resistance and beam patterns for the case N=4. In this thesis we present graphs for the 
case N=6. It will be shown that the six resonator case is superior to the four resonator 
case, and is in fact the best choice. In the following sections we present results for the 


planar and the conical configurations. 
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A. THE PLANAR CONFIGURATION 


1. THE GEOMETRY 
Figure 3-1 shows the geometry of the planar circular configuration for the 
case of six passive acoustic resonators. The radius of the circle is R and the transducer 
is in the plane of the resonators (x-y plane). We are interested in dimensionless 
distances, so the radius of the circle will be referred to in terms of wavenumber times 
the radius, KR. The dimensionless radius of each resonator is ka and the dimensionless 


radius of the transducer is kay . We assume a value for the transducer radius of 





Figure 3-1 Geometry for the planar configuration with six resonators. The plot 1s drawn 
to scale for resonators with ka = 0.1. 
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kay=0.2. Later we will demonstrate that the results are independent of the value chosen 


over the range of interest. 


2. GAIN IN RADIATION RESISTANCE 

Figure 3-2 shows a three-dimensional surface plot and contour plot of the 
gain in radiation resistance versus ka and kR for six air bubbles, at 1 atm ambient 
pressure, in the configuration of Figure 3-1. The gain has a peak value of 1.95 at ka = 
0.0138 and KR = 2.87. The increasing value of gain as КК gets small is not interesting 
because of the expected divergent behavior [Ref. 4]. Inspection of the gain plot reveals 
that the peak in gain 1s very narrow in the ka-direction. The reason is the very sharp 
peak in the tuning curve for a resonant bubble. The Q of a resonant bubble is about 
ee 

Тһе gain in radiation resistance was also plotted for a bubble-like passive 
acoustic resonator whose resonant value of ka is 0.1. This is shown in Figure 3-3. An 
air bubble at 50 atm hydrostatic pressure will have ka,=0.1. It is also expected that a 
more depth-independent resonator such as a thin-shell could be constructed with ka,=0.1. 
The Q for this case is about 10 (as compared to about 70 for an air bubble at 1 atm 
ambient pressure.) Figure 3-3 shows a much broader range of ka values with substantial 
(z1.5) gain. The implication of these observations is that using air bubbles at 1 atm 
ambient pressure would not be practical, but that using resonators with ka,=0.1 could 
result in a practical application. 

Table 3-1 shows numerical results for the planar configuration for two 


through ten resonators at 1 atm ambient pressure. The calculations were performed using 
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Figure 3-2 Gain in radiation resistance versus ka and KR for the planar configuration 
with six air bubbles at 1 atm ambient pressure (a) 3D surface plot (b) contour plot 
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Figure 3-3 Gain in radiation resistance versus ka and KR for the planar configuration for 
six resonators with ka, = 0.1 (a) 3D surface plot (b) contour plot 
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Table 3-1 Planar configuration results for air bubbles at 1 atm ambient pressure. 


ашат 
90.0 
[os 1.95490 90.0 
0.05 aaa] 
[0975 | 010138355 | 2.86790 | 1.95482 | 0.47270 | 900 _ 
LAN EMT EE 
0.0138355 | 2.86927 | 1.95350 | 0.464611 
КТЕТІТІ-ЕЕЕГІТЕ ЕСІТЕ ШІСІ БЕСТЕН 
_ол5 | 0.0138342 | 292865 | 1.89464 | 0367996 | 900 


Table 3-2 Planar configuration results for six resonators allowing Ka, to vary. 
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Mathematica with a precision of about ten significant digits (the working precision was 
twenty significant digits). Fewer digits are shown in the table because they are not of 
general interest. Note that the table indicates that there is no additional increase in gain 
beyond six resonators. This is in part the basis for our claim that, for the circular 
configurations, six resonators is the best choice. Table 3-1 also reveals that when the 
gain reaches its maximum value, the total volume velocity of the resonators (N. | U” |) 
remains constant with increasing numbers of resonators. The phase of the volume 
velocity of each resonator (/U"”) is 90° relative to the transducer. We will make use 
of this fact in the conical configuration. 

Recall that the calculations for Figures 3-1 and 3-2 and Table 3-1 are for 


kay=0.2. Table 3-2 shows the results for six resonators where we allow kay to vary. 


_2 | 00096243 | 265 | 126 | 0.388 | 076 |900 |12. 
13 | 0.093435 | 2.70 | 159 | 0,576 | 173 |900 |42. 
4 | 0.997272 | 2.82 | 1.87 | 0.639 | 2% |900 |62 
5 | 0.100877 | 2.86 | 1.94 | 0556 | 278 |900 |67. 


_6 | owo117 | 2.87 | 1.95 | 0467 | 209 |90 67 
7 | o.101548 | 2,87 | 95 | 0400 | 280 | 90.0 | 67 
св | 0101813 | 287 | 195 | 0350 | 280 | 900167. 
_9 | 101949 | 2.87 | 1.95 | 0.311 | 20 |900 |67 


10 | 0.101981 | 2.87 | 1.95 | 0.280 2.80 | 90.0 


Table 3-3 Planar configuration results for resonators with ka,=0.1. 
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The results show that, over a large range of values, the gain in radiation resistance is 
independent of the value of ka, chosen. The numerical results for resonators with 
ka,=0.1 are shown in Table 3-3. The values of KR, Gain, U“, and directivity index (DD) 
remain the same as in Table 3-1. The values of Xa are different but it can be shown that 
the ratio ka/ka, is the same. This indicates that over a large range of values, the gain in 


radiation resistance is independent of the value of ka, for bubble-like resonators. 


3. BEAM PATTERNS AND DIRECTIVITY 
The three-dimensional beam pattern can be found by plotting equation 2-24 


over all spherical angles. Figure 3-4 1s the three-dimensional far-field beam pattern for 
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Figure 3-4 Three-dimensional beam pattern for planar configuration with six 
resonators at 1 atm ambient pressure. 
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the planar configuration with six air bubbles at 1 atm ambient pressure. The three- 
dimensional (3D) beam pattern provides the reader with an idea of what the pressure 
field will look like in 3D space without having to resort to multiple 2D plots in various 
planes. Now that we know what the beam pattern looks like we can better utilize our 
polar plots to analyze the results. From the 3D plot alone, we can see that there is some 
significant directivity along the acoustic axis (in this case the positive z-axis). We can 
also see that the pressure field is symmetric about the z-axis and the x-y plane. 

The usual method of presenting a beam pattern is the dB polar plot. Figure 


3-5 is a polar plot of the beam pattern for the planar configuration with six resonators. 


aa as 





Figure 3-5 Polar plot of the beam pattern for the planar configuration with six 
resonators at | atm ambient pressure. 
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The zero dB reference is the pressure that would exist due to the transducer if no passive 
elements were present. The dashed line represents the sound pressure level that would 
exist due to a spherical source radiating the same acoustic power as the transducer does 
in the presence of the passive elements. The directivity index, defined by equation 2-26, 
can be read directly from Figure 3-5 as the difference in dB between the sound pressure 
level of the transducer in the presence of passive acoustic radiators and that of the 
spherical source radiating the same acoustic power. The directivity index in this case is 
6.7 dB. Tables 3-1 and 3-3 list the directivity index for the planar configuration for two 


through ten resonators at 1 atm ambient pressure and for ka,=0.1. 


В. THE CONICAL CONFIGURATION 


1. THE GEOMETRY 

The relative volume velocity of the passive acoustic resonators was noted in 
the last section to be ninety degrees out of phase with the transducer volume velocity. 
We also noted that the far-field beam pattern was symmetric about the x-y plane. For 
a directional source, we usually want the sound energy to be directed along a single 
acoustic axis with any side or back lobes minimized. We make use of the phase of the 
resonators’ relative volume velocity in the conical configuration to improve the directivity 
of the system. In the conical configuration the transducer is displaced one-quarter 
wavelength out of the plane of the resonators. Figure 3-6 shows the geometry of the 


conical configuration for the case of six passive acoustic resonators. The dimensionless 


2 





Figure 3-6 Geometry of the conical configuration with six resonators. 


radius of the circle of resonators is KR and the dimensionless height of the transducer is 
КН = 1/2. 
2. GAIN IN RADIATION RESISTANCE 

Surface and contour plots of the gain in radiation resistance for the conical 
configuration with six resonators at 1 atm ambient pressure are shown in Figure 3-7. 
The plots are similar to the planar configuration plots except that the gain is not as large. 
The peak in gain seen in Figure 3-7 is 1.82 at ka 20.0138 and KR—2.69. The value of 
ka is 1.005 times the resonant value for a single bubble at 1 atm ambient pressure and 
is approximately the same as the value of ka for maximum gain in the planar 


configuration. As in the planar configuration, the peak in gain is very narrow in the 
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Figure 3-7 Gain in radiation resistance versus ka and kR for the conical configuration 
with six resonators at 1 atm ambient pressure (a) 3D surface plot, (b) contour plot 
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Figure 3-8 Gain in radiation resistance versus ka and KR for the conical configuration 


with six resonators for ka, = 0.1 (a) 3D surface plot (b) contour plot 


32 


ka-direction. Figure 3-8 shows the surface and contour plots for the conical 
configuration where the resonators are assumed to have a value of ka,=0.1. Once again, 
the peak has broader support in the ka-direction for ka,=0.1. 

The numerical results for the conical configuration are given in Tables 3-4, 


3-5, and 3-6. Table 3-4 gives the results for air bubbles at 1 atm ambient pressure for 
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Table 3-4 Conical configuration results for air bubbles at 1 atm ambient pressure. 


two through ten resonators. The conical configuration, like the planar configuration, 
does not provide any additional gain in radiation resistance for more than six resonators. 
Table 3-5 gives the results for six resonators in a conical configuration where we allow 
ka, to vary from the usual value of 0.2. The gain in radiation resistance does not vary 


significantly with the value of ka, chosen. Table 3-6 presents the results for the conical 
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Table 3-5 Conical configuration results for six resonators allowing kay to vary. 
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Table 3-6  Conical configuration results for passive acoustic resonators with 
ka,=0.1. 
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configuration with resonators having a value of ka,=0.1. The values of kR, Gain, U””, 
апа D/ remain about the same as in Table 3-4. While the values of ka differ, the ratio 
ka/ka, is the same for both the 1 atm bubble and the ka,=0.1 resonator. 

Note that, for four or more resonators, the planar configuration has a larger 
gain in radiation resistance than the conical configuration with the same number of 
resonators. Calculations show that for four or more resonators, the maximum gain 


always occurs when the transducer is in the plane of the resonators. 


3. BEAM PATTERNS AND DIRECTIVITY 

Figure 3-9 shows the three-dimensional far-field beam pattern, defined by 
equation 2-24 and plotted over all spherical angles. The improvement over the planar 
configuration is obvious. We traded a small fraction of the gain in radiation resistance 
for a large increase in directivity. In addition to concentrating sound energy along an 
acoustic axis, the conical configuration also has srnaller side and back lobes than the 
planar configuration. Figure 3-10 is the polar dB plot of the conical configuration with 
six resonators at 1 atm ambient pressure. From Figure 3-10 we note that the conical 
configuration concentrates 3 dB more sound pressure along the acoustic axis than the 
planar configuration. This represents a doubling of the acoustic intensity. 

The directivity index 1s determined from Figure 3-10 as a described in the last 
section. The directivity index for six resonators is 8.4 dB. Tables 3-4 and 3-6 list the 
directivity index for two through ten resonators at 1 atm ambient pressure and for 
ka,=0.1. The directivity index, like the gain, does not increase with more resonators 


past its value for six resonators. The high gain in radiation resistance and high 
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Figure 3-9 Three-dimensional beam pattern for the conical configuration with six 
resonators at | atm ambient pressure. 


directivity index make the conical configuration the best choice for a system design. As 
we will see in the following chapters, the 3D beam pattern for the conical configuration 


is the most desirable of the configurations considered. 
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Figure 3-10 Polar plot of the beam pattern on a dB scale for the conical 
configuration with six resonators at 1 atm ambient pressure. 
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IV. LINE CONFIGURATION TYPE I 


In this chapter and the next two, we extend the work done by Ellsworth by 
considering linear configurations of passive acoustic resonators. The motivation is to 
investigate whether a linear configuration offers any improvement in power radiated and 
directivity compared to a circular configuration. We present results for six resonators 
for easy comparison to the results in the last chapter. The first line configuration 
considered, termed line configuration type I, consists of a line of equally spaced 
resonators on only one side of the transducer. In the next chapter, we consider an equal 
number of resonators on either side of a transducer in a line (line type II). In Chapter 


V], we will consider a Yagi-Uda [Ref. 5] array of resonators (line type III). 


A. THE GEOMETRY 

Figure 4-1 shows the line configuration type I with six resonators. The transducer 
is at the origin of the coordinate system and the resonators are along the positive x-axis. 
The dimensionless distance from the transducer to the first resonator is AL. The 
dimensionless distance between resonators is k/. The transducer radius is kag=0.2. Each 
passive acoustic resonator has radius ka,, where the subscript 7 indicates the individual 


resonator. 
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Figure 4-1 Geometry of the line configuration type I with six resonators. 
B. GAIN IN RADIATION RESISTANCE 


Figure 4-2 shows the gain in radiation resistance presented to the transducer in the 
presence of a system of six paSsive acoustic radiators arranged in line configuration type 
I. The gain in radiation resistance is plotted versus kL and kl with ka held constant at 
the value that provides maximum gain for a finite value of AL. The maximum considered 
is the one at AL=3.6 and k/=2.5, for which ka=0.01373. The gain at this point is 1.22 
(compared to the conical configuration with a gain of 1.82 for six resonators). The peaks 
that occur in the gain as KL goes to zero were only briefly considered. This 15 discussed 
in the next section. A plot of the gain in radiation resistance versus AL and kl for 
ka=ka,=0.1 1s indistinguishable from Figure 4-2 and is not presented here. We expect 


that, as for the circular configurations, the ka, —0.1 case 1s preferable because the range 
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Figure 4-2 Gain in radiation resistance versus kL and kl for the line configuration 
type I with six resonators at 1 atm ambient pressure (a) 3D surface plot (b) contour plot 


of acceptable values of ka, i.e. which yield a gain greater than a certain value, would 
be greater. 

Tables 4-1 and 4-2 present results for the line configuration type I for one through 
fourteen resonators at 1 atm ambient pressure and with ka,=0.1 respectively. Notice that 
the gain remains fairly constant at about 1.2. The relative volume velocities listed are 
those for the resonator closest to the transducer. As before, the primary difference 
between the two tables is the values of ka. It can be shown that the ratio ka/ka, is the 


same for both. 


C. BEAM PATTERNS AND DIRECTIVITY 

Figure 4-3 shows the three-dimensional far-field beam pattern for the line type I 
configuration with six resonators at 1 atm ambient pressure. The acoustic axis is along 
the negative x-axis. Figure 4-4 shows a polar plot of the beam pattern. The directivity 
index, found from Figure 4-4 or Table 4-1, is 4.6. The line configuration type I can be 
used to increase the power output and directivity of an underwater transducer; however, 
the conical configuration results are clearly better. 

In the last section, we described divergent peaks in the plot of gain in radiation 
resistance for vanishing AL. For one of these peaks, for example, the gain at AL=1 and 
kl=2.33 is 1.82. We plotted a beam pattern using these values in order to better 
understand these peaks in the gain plot. The directivity index was found to be DJ=11.0. 
Figure 4-5 is a polar plot of the resulting beam pattern. The acoustic axis is now the 


positive x-axis. This is an indication that the peaks as kL goes to zero correspond to a 
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Table 4-1 Line configuration type I results for air bubbles at 1 atm ambient 


pressure. 
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Table 4-2 Line configuration type I results for resonators with ka,=0.1. 
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Figure 4-3 Three-dimensional beam pattern for the line configuration type I with 
six resonators at 1 atm ambient pressure. 


condition on the values of kl that provide constructive interference. Itis probably related 
to the total length of the aperture. Further investigation is required to fully understand 
the occurence of these gain peaks and to determine whether a practical system can utilize 
this phenomenon. If exploitable, this system provides the best directivity of any analyzed 


in this research. 
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Figure 4-4 Polar plot of the beam pattern for the line configuration type I with six 
resonators at 1 atm ambient pressure 





Figure 4-5 Polar plot of the beam pattern for the line configuration type I with six 
resonators at ] atm ambient pressure for the divergent peak. AL=1, kl=2.33 
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V. LINE CONFIGURATION TYPE II 


A. THE GEOMETRY 

Figure 5-1 shows the geometry for the line configuration type II with six 
resonators. In this configuration, half of the resonators are on one side of the transducer 
and half are on the other side. The dimensionless distance from the transducer to thefirst 
resonator on either side is AL. The dimensionless inter-resonator spacing is kl. The 


transducer is assumed to have a dimensionless radius of ka,=0.2. 





Figure 5-1 Geometry of the line configuration type II with six passive acoustic 
resonators. 
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В. GAIN IN RADIATION RESISTANCE 

The gain in radiation resistance versus AL and kl for a fixed ka is shown in Figure 
5-2. The value chosen for ka is the value that provides maximum gain at the local 
maxima at KL=3.0 and kl=2.2, for which ka=0.01376. The gain at this point is 1.56. 
As with the line type I, the gain plot for line type II diverges as KL goes to zero. We 
have not investigated this behavior in detail. Tables 5-1 and 5-2 give the numerical 
results for line type II for two through twenty air bubbles at 1 atm ambient pressure and 


for resonators with ka,=0.1 respectively. 
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Table 5-1 Line configuration type II results for air bubbles at ] atm ambient 
pressure. 
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Figure 5-2 Gain in radiation resistance versus kL and E for the line configuration 
type II with six resonators at 1 atm ambient pressure (a) 3D surface plot (b) contour plot 
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As seen in Table 5-1, the gain in radiation resistance does not reach a constant value 
(over the range of N considered). For twenty resonators the gain is 2.17. This is the 
highest gain calculated in this research. The values of KI and kL are the same in Tables 
5-1 and 5-2. The difference is in the values of ka. It can be shown that the ratio ka/ka, 


is the same in both cases. 








(y| e [uu | Goin 1010) р 

з [кше за [ae [aa Tn [or 

4 | 0.100526 | 1.89 | 2.98 | 1.42 | 0.633 | 862 | 29 — 

[6 | 00996309 | 2.24 | 2.99 | 1.56 | 0.756 | 809 | 41 

_8 | 0.0990346 | 2.38 | 3.06 | 1.68 | 0857 | ssa | so 
18 | 00974035 2.70 | 3. 17 DELI 1.08 Sad 


E NO 


Table 5-2 Line configuration type II results for resonator with ka,=0.1. 











C. BEAM PATTERNS AND DIRECTIVITY 
Figure 5-3 shows the three-dimensional far-field beam pattern for the line 
configuration type II with six resonators at 1 atm ambient pressure. The beam pattern 


for this configuration is not particularily interesting. The acoustic axis is along the 
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Figure 5-3 Three-dimensional beam pattern for the line configuration type II with 
six resonators at 1 atm ambient pressure. 


positive and negative axes and there is some directivity. Figure 5-4 shows the dB beam 
pattern. The directivity index (D/), from Figure 5-4 or Table 5-1, is 4.9 dB. Unlike the 
circular configurations, the DI does not reach a constant valve (over the range of N 
considered). For twenty resonators, the line type II provides a directivity index of 6.5 
dB. 

Figure 5-5 shows the beam pattern for this configuration when using the value of k/ 
corresponding to the peak at kL=1. The value of k/ is 1.72, the directivity index is 8.0 


dB, and the gain is 2.15. The large gain and high directivity makes this system 
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Figure 5-4 Polar plot of the beam pattern for the line configuration type II with six 
resonators at 1 atm ambient pressure 





Figure 5-5 Polar plot of the beam pattern for the line configuration type II with six 
resonators at 1 atm ambient pressure for the divergent peak. KL-1, k/=1.72 
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potentially useful for practical application. However, further investigation is required 
to understand this behavior and to determine whether a practical system can utilize the 
large gain in radiation resistance at small values of KL. We provide the result as a basis 


for this investigation. 
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VI. AN ACOUSTIC YAGI-UDA ARRAY 


The Yagi-Uda electomagnetic antenna system [Ref. 5] uses parasitic antenna 
elements around a feed antenna element in order to direct a travelling wave. It depends 
on mutual coupling between the elements. The analogy to the underwater transducer 
with passive acoustic resonators is direct. The transducer is the feed element and the 
passive acoustic resonators are the parasitic elements. In the underwater case, the 
elements are coupled by the sound field. We are basically trying to do the same thing 
that antenna farm designers have been doing for years with radio waves. Once you get 
past the obvious analogy, the similarities end. The acoustic Yagi-Uda array described 
here is more closely related to the work described earlier in this paper than to the 
electromagnetic antenna system. We will alternately refer to the Yagi-Uda array as line 
configuration type III. In the sections that follow, we will develop the theory and then 


present numerical and graphical results. 


A. THEORY 

In the design of the Yagi-Uda configuration (line type III), we are interested in 
directing the sound energy by choosing resonators that will produce constructive 
interference in one direction and destructive interference in the opposite direction. We 


also impose the condition that the gain in radiation resistance be not less than one. 
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For constructive interference, we require that the phase of the pressure wave 
outgoing from a resonator equal the phase of the pressure wave incident at the resonator 
from the transducer. This constraint can be met only if the resonator is stiffness 
controlled. We then place an additional constraint on the system requiring that the gain 
be one or greater. For this condition, we must have the phase of the pressure wave 
incident at the transducer from a resonator equal to the phase of the pressure wave 
outgoing from the transducer. These two constraints require (1) the resonators must be 


stiffness-controlled, and (2) the dimensionless distance from the transducer must satisfy 
kd = mm, т=|,2,%,... (6-1) 
where d, = the distance from each resonator to the transducer 


and k wavenumber in the medium. 


For destructive interference, we require that the phase of the pressure wave 
outgoing from a resonator equal «x plus the phase of the pressure wave incident at the 
resonator from the transducer. This constraint can only be met if the resonator is mass- 
controlled. We then use the constraint that the gain be one or greater. These two 
contraints require (1) the resonators must be mass-controlled, and.(2) the dimensionless 


distance from the transducer must satisfy 


kd = Int oe nepos wm (6-2) 


aft 
b 2 


We force the resonators to be approximately stiffness-controlled yet still close to 


resonance by setting ka =ka, - 1/2(ka,)’. We force the resonators to be approximately 
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mass-controlled and likewise close to resonance by setting ka = ka, + 1/2(ka,)°. These 
values result in the resonators being driven at a frequency which is one resonant half- 
width above and below their resonant frequency (the resonant half-width is 1/2(ka,)’). 


Figure 6-1 is a diagram of the line type III configuration. We show the distances 
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Figure 6-1 Diagram of the Yagi-Uda configuration (line type III) for six resonators. 


in terms of wavelengths and indicate which resonators are mass- and stiffness- controlled. 


In the next section we present the results of this research. 


B. RESULTS 

Computer simulations of the line type III configuration indicate that the analysis of 
the theory in the previous section is correct. We did not perform any optimization of 
parameters on this configuration as we did for previous configurations. The results 
presented here are for the distances and radii described in the previous section. 

Figure 6-2 shows a three-dimensional far-field beam pattern for a line type III 
configuration of six resonators. Тһе acoustic axis is along the positive x-axis, as 


expected. Figure 6-3 shows the polar plot of the beam pattern. Table 6-1 lists the gain 
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Table 6-1 Results for the Yagi-Uda configuration (line type III) 

in radiation resistance and directivity index for two through ten resonators. The gain for 
six or more resonators is about one and the directivity index is about 6.5 dB. Figure 6-4 
is the beam pattern and Figure 6-5 the polar plot for the case of ten passive acoustic 
radiators, and is provided for comparison. As indicated in Table 6-1, there is no major 
difference in either the gain in radiz.tion resistance or the directivity index for the N=6 


and N=10 cases. 


D» 





Figure 6-2 Three-dimensional beam pattern for line type III configuration with six resonators 





Figure 6-3 Polar plot of the beam pattern for the line type III configuration with six resonators 
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Figure 6-4 Three-dimensional beam pattern for line type III configuration with ten resonators 





Figure 6-5 Polar plot of the beam pattern for the line type III configuration with ten resonators 
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VII. DISCUSSION 


We have analyzed five configurations of resonators and provided many graphs and 
tables of results. Table 7-1 is a comparison of the results for all configurations. The 
measures of effectiveness that we chose for comparison are (1) gain in radiation 


resistance, (2) directivity index, and (3) the far-field beam patterns. Table 7-1 lists the 






Planar | Conical 

цат = al 
А е ж. "n E. 
I EO JE "Toe" ІН Kë D 1.26 | 

120 E 
: 
[s|195| 67 | 132 |84| 122 | 52 | 168 | «9 | 097 Les 
[i| 196] 67 | 182 | 84 | 100 | 54 | 178 | 55 | 099 Les 


Table 7-1 Comparison of the results for all configurations 












data for the first two measures of effectiveness and Figure 7-1 provides a comparison of 
the far-field beam patterns on a single plot. Recall that O dB represents the transducer 
with no resonators present. We considered the case of six resonators throughout the 
investigation as a basis for comparison of the three-dimensional beam patterns. The 
N =6 case is highlighted in Table 7-1. The planar configuration has the best gain, almost 


double, and the conical configuration has the best directivity index, 8.4 dB. Figure 7-1 
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Acoustic Axis 





=т= = Line Ill 


Figure 7-1 Polar plot of the far-field beam pattern for all configurations for six 
resonators. 
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shows that the conical configuration also has the best front-to-back discrimination along 
the acoustic axis. Of the results presented, clearly the conical configuration is the best 
choice for a real system. The linear configurations are generally worse than the circular 
configurations, but the line type III configuration is of particular interest. The addition 
of six passive acoustic resonators in a Yagi-Uda antenna-like configuration greatly 
improves the directivity of the transducer. However, since the passive resonators are not 
driven at their resonant frequency, this configuration does not provide any power gain. 

Further investigation of the divergent peaks in gain for the line configurations type 
I and II is warranted. The gain diverges as the distance, L, to the first resonator tends 
toward zero. It may be possible to exploit regions of finite KL for which this behavior 
occurs. Certain values of resonator spacing satisfy a resonant condition related to the 
total length of the aperture which apparently results in substantial gain. We presented 
beam patterns for a system operating on one of these divergent gain peaks and found that 
the directivity and gain of the system were better than any other system analyzed. 
Further investigation into this behavior 1s encouraged. 

We have shown that resonators with a value of ka, — 0.1 provide identical results 
in terms of gain and directivity as do bubbles at 1 atm ambient pressure. We also 
demonstrated that these resonators are not as sensitive to the value of ka, that is, 
substantial gain (21.5) may be achieved over a reasonable range of operating 
frequencies. This is understood as follows. For an air-filled bubble at 1 atm ambient 
pressure, Q=1/ka, is about 70. For a system to be feasible, the dimensions of each 


resonator must be matched to substantially better than one part in Q. It is probably not 


60 


reasonable to expect that a collection of bubbles can be matched closely enough in radius 
for a system at 1 atm ambient pressure to be practical. It may be possible to build, 
however, a system of bubble-like resonators such as thin-walled shells that will have a 
ka, of about 0.1 and therefore a Q of about 10. 

Table 7-2 describes the physical size of a system for various configurations of six 


resonators as described in this thesis. The overall dimensions of the system are the same 


| 


Radius of (cm): 
Resonator 33 240 3.3 24 0.33 


Aperture size of (m): 
Planar 137 Ise 137 


128 


Line Type I 9913 3.83 


Line Type I EL v s» 
Line Tye | as | as | an 





Table 7-2 Physical size of the system of resonators for 10 Hz, 100 Hz, and 1 kHz for 
various configurations of six resonators. 


for bubbles at 1 atm ambient pressure and for ka,=0.1 resonators. The size of the 
resonator increases for larger values of ka,. The conical configuration, chosen as the 
best in terms gain and directivity, has the smallest physical dimensions. The aperture 


size (diameter for circular and total length for linear configurations) for the conical 
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configuration at 100Hz is about 13 


have a radius of 24 cm. 


VIII. CONCLUSION 


А. SUMMARY 

We developed the equations for the network analysis of a transducer in the presence 
of passive acoustic radiators. These equations were implemented in computer programs 
and applied to five configurations of resonators. We first considered the two 
configurations previously examined by Ellsworth. We provided results for the case of 
six resonators and developed three-dimensional beam patterns. We then examined two 
linear configurations and provided results. Next, we developed the idea of a Yagi-Uda 
antenna-like array of passive acoustic resonators. Finally, we conducted a comparison 


of all five configurations. 


B. CONCLUSIONS 

We found that the conical configuration of six resonators is the best choice in terms 
of the gain in radiated power output, directivity along an acoustic axis, and minimum 
required number of resonators. A linear configuration could be utilized to improve the 
gain and directivity of a transducer, but not as well as the conical configuration. We 
found that the acoustic Yagi-Uda antenna configuration can be used to improve the 


directivity of an underwater transducer. 
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C. RECOMMENDATIONS 

The results of the research presented in this paper are promising. We recommend 
that further work be done to (1) design a resonator with ka, = 0.1, (2) build a system 
of resonators and verify the theory and results in this paper, (3) study the Yagi-Uda 
configuration in more detail including optimizations to improve the gain , and (4) study 
the divergent peaks in gain as kL goes to zero for the linear configurations to determine 
whether a practical system can be designed to utilize the higher theoretical gain and 


directivity. 


APPENDIX 


The following Mathematica program code is an example of how programs were 
written to perform the calculations in this paper. The example consists of the standard 
definitions for the conical configuration and code to create plots as found in this 
document and calculate the information found in tables for the conical configuration. 


The lines in bold face type are comments. 


Distance, Radius, and Area 


ClearAll[ka,S,kl,krange] 


Ка[0] : = ka 
ka[n_] := kai 
S[n_] := ka[n] * 2 
kl[n_,n ]:=0 


kl[n_,0] := Sqri[kR*2 + kH”2] 

kI[0,n_] := Sqrt[kR*2 + kH”2] 

kl[n_,m_] := КК Sqrt[2-2*Cos[(2 Pi/Num) (n-m)]] 

krange[0,the_,ph_] := kr 

krange[n ,the ,ph ]:- Sqrt[ 
(kr Sin[the] Cos[ph] - kR Cos[(2 Pi n)/Num)*2 + 
(kr Sin[the] Sin[ph] - kR Sin[Q Pi n)/Num))"2 + 
(kr Cos[the] - kH)”2] 


Impedances 


ClearAll[Za,Zar,Dfc] 

Dfc[n_] := Exp[j ka[n]] / (1+ka[n)) 

Za[n_} := - I (S[O]/S[n]) ((3 gam Po)/(p c^2 ka[n])) 

Zar[n_,n_] := (S[O]/S[n]) (1 ka[n]/(1 +1 ka[n])) 

Zar[n_,m_] := Dfc[n] Zar[m,m] (ka[m]/kl[n,m]) Exp[ 
I (ka[m] - ki[n,m])] 
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Gain in radiation resistance 


ClearAll[ZB,ZC, Gain] 

ZB[n_] := Zar[n,0] 

ZC[n_,n_] := Zar[n,n] + Za[n] 

ZC[n ,m ] :— Zar[n,m] 

Gain :— Re[Zar[0,0] - Array[ZB, Num] . Inverse[Array[ 
ZC,{Num,Num}]] . Array[ZB,Num]]/Re[Zar[0,0]] 


Relative velocities and field pressure 


Clear[theta,phi,the,ph,n] 
ClearAll[BP,dBBP,Prad, Urel] 
Urel := - Inverse[Array[ZC,íNum,Numj]] . Array[ZB,Num]; 
Prad[0,the_,ph_] := Zar[0,0] (ka[O]/krange[O,the,ph]) Exp[ 
I (ka[0]-krange[O, the, ph])] 
Prad[n ,the ,ph ] :- Zar[n,n] Urel[[n]] (ka[n]/ 
krange[n,the,ph]) Exp[I (ka[n]-krange[n,the,ph])] 
BP[the ,ph ] :- Abs[Sum[Prad[n,the,ph], (n,0, Numj]]/ 
Abs[Prad[0, the, ph]] 
dBBP[the_,ph_] := 20 Log[10,BP[the,ph]] 


Constants 


Clear{kH,kai,kaO,kR,kr,kar,gam,Po,c,p] 
gam=1.402; 

Po=1.013 10°S; 

c= 1500; 

p=998; 

kar = Sqrt[(3 gam Ро)/(р с^2)] 

Num =6 

КН =P1/2 

Ка0-0.2 


Create surface plot of Gain versus ka and kR 
plot! =Plot3D[Gain, {kR,1,5}, {kai,0.9 kar,1.1 kar}, 


PlotPoints- > 40,PlotRange- > {0,2}, 
AxesLabel- > ("kR","ka","Gain")] 
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Create contour plot 


plot2 2 Show[ContourGraphics[plot1], 
Frame- » True, 
FrameLabel- » ("kR","ka"!, 
RotateLabel- » True, 
ContourSmoothing- > 3, 
ContourShading- > False, 
AspectRatio- > 1, 
Contours- > {0.2,0.4,0.6,0.8,1.0,1.2,1.4,1.6,1.83] 


Find Maximum Gain 


Clear[kai, kR] 
{maxgain, {kaiRule,kRrule}}= 
FindMinimum[-Gain, (kai, {kar, 1.002 kar}}, 
ERO TON 
kai — kai/. kaiRule 
kR=kR/.kRrule 
MaxGain =-maxgain 


Find relative volume velocities 
N[(Abs[Urel[[1]]], Arg[Urel[[1]]] 180/Piy] 
Find directivity index 


thAxis = 0; (* theta on acoustic axis *) 
phAxis = 0; (* phi on acoustic axis *) 
DI = dBBP[thAxis,phAxis] - 

20 Log[10,Sqrt[MaxGain]] //N 


Create 3D beam pattern 


BP2[th_,ph_]=Evaluate[BP[th, ph]]; 

Needs["Graphics'ParametricPlot3D*"] 

plot3 = ParametricPlot3 D[BP2[theta, phi] 
(Sin[theta] Cos[phi], 
Sin[theta] Sin[phi], 
Cos[theta]}, 
(theta,0, Pi, Pi/40}, 
phi,0, 2 P1, P1/20]] 
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Create 2D (dB) beam pattern in z-x plane 


data=Table{N[dBBP[theta,0]], 
{theta,0,2 Pi, Pi/40}); 
mystuff=Table[(data[{i]} + 12) 
(N[Sin[G-1) Pi/40]], N(Cos((-1) Pi/40]]]), 
ІРІ, ӨР 
myplot=ListPlot[mystuff, PlotJoined- > True, 
AspectRatio- > Automatic] 
The following requires a separately created template and label: 
Show[template,label,myplot, 
AspectRatio- > Automatic, 
Axes- > False, PlotRange- > All] 
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